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Superconductivity and magnetism in the non-centrosymmetric heavy fermion compound 
CePtsSi and related materials are theoretically investigated. On the basis of the random phase 
approximation (RPA) analysis of the extended Hubbard model, we describe the helical spin fluc- 
tuation induced by the Rashba-type anti-symmetric spin-orbit coupling and identify two stable 
superconducting phases with either the dominant p-wave (s+P-wave) symmetry or the d-wave 
(p+D+/-wave) symmetry. The effect of the coexistent antiferromagnetic order is investigated in 
both states. The superconducting order parameter, quasiparticle density of state, NMR 1/TiT, 
specific heat, anisotropy of H C 2, and possible multiple phase transitions are discussed in detail. 
A comparison with experimental results indicates that the s+P-wave superconducting state is 
likely realized in CePtsSi. 
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1. Introduction 

The discovery of superconductivity in materials with- 
out an inversion center 1 ' 2 has initiated intensive research 
on a new aspect of unconventional superconductivity. 
Several new non-centrosymmetric superconductors with 
unique properties have been identified among heavy 
fermion systems such as CePt 3 Si, 1,2 Ulr, 3 CcRhSi 3 , 4 ' 5 
CeIrSi 3 , 6 ' 7 CeCoGe 3 8 and others like LiaPd^Ptj^B, 9 
Y 2 C 3 , 10 Rh 2 Ga 9 , Ir 2 Ga 9 , n Mg 10 Ir 19 B 16 , 12 > 13 Re 3 W, 14 
boron-doped SiC 15 and some organic materials. 16 

One immediate consequence of non- 
centrosymmetricity is the necessity for a revised 
classification scheme of Cooper pairing states, as parity 
is not available as a distinguishing symmetry. Super- 
conducting (SC) states are considered as a mixture of 
pairing states with different parities or, equivalently, 
the spin configuration is composed of both a singlet 
component and a triplet component. The mixing of 
spin singlet and spin triplet pairings is induced by 
the anti-symmetric spin-orbit coupling (ASOC). 17 
Recent theoretical studies led to the discussion of 
interesting properties of a non-centrosymmetric su- 
perconductor, such as the magnetelectric effect, 18-21 
anisotropic spin susceptibility 17 ' 21-29 accompanied by 
the anomalous paramagnetic depairing effect, 29 anoma- 
lous coherence factor in NMR 1/TiT, 20 ' 30 anisotropic 
SC gap, 28 ' 30-33 helical SC phase, 29 ' 34-39 Fuldc-Ferrel- 
Larkin-Ovchinnikov (FFLO) state at zero magnetic 
field, 40 various impurity effects, 41-44 vortex state, 45 ' 46 
and tunneling/ Josephson effect. 47-52 

Non-centrosymmetric heavy fermion superconductors, 
i.e., CePt 3 Si, Ulr, CcRhSi 3 , CeIrSi 3 , and CeCoGe 3 , are 
of particular interests because non-s-wave superconduc- 
tivity is realized owing to strong electron correlation 
effects and magnetism has an important effect on the 
superconducting phase. However, the relation between 



magnetism and superconductivity has not been theoreti- 
cally studied so far, except in studies refs. 28, 53, and 54. 
Here, we extend our previous study 28 and investigate the 
pairing state arising from magnetic fluctuation in detail. 

Another aim of this study is to elucidate the effects of 
antiferromagnetic (AFM) order on the SC phase. Inter- 
estingly, all presently known non-centrosymmetric heavy 
fermion superconductors coexist with magnetism. We 
have shown that some unique properties of CePt 3 Si at 
ambient pressure can be induced by the AFM order. 28,29 
In this study, we analyze this issue in more detail. 

Among non-centrosymmetric heavy fermion supercon- 
ductors, CePt 3 Si has been investigated most extensively 
because its superconductivity occurs at ambient pres- 
sure; others supcrconduct only under substantial pres- 
sure. Therefore, we focus here on CePt 3 Si. We believe 
that some of our results are qualitatively valid for other 
compounds too. In CePt 3 Si, superconductivity with T c <~ 
0.5K appears in the AFM state with a Nccl tempera- 
ture Xn = 2.2K. 1 The AFM order microscopically coex- 
ists with superconductivity. 55 ' 56 Neutron scattering mea- 
surements characterize the AFM order with an ordering 
wave vector Q = (0, 0, it) and magnetic moments in the 
ab- plane of a tetragonal crystal lattice. 57 The AFM or- 
der is suppressed by pressure and vanishes at a critical 
pressure P c ~ 0.6GPa. Superconductivity is more ro- 
bust against pressure and therefore a purely SC phase is 
present above the critical pressure P > 0.6GPa. 58-60 

The nature of the SC phase has been clarified by 
several experiments. The low-temperature properties 
of thermal conductivity, 61 supcrfluid density, 62 specific 
heat, 63 and NMR 1 /7\T 64 indicate line nodes in the gap. 
The upper critical field H C 2 ~ 3 — 4T exceeds the stan- 
dard paramagnetic limit, 1 which seems to be consistent 
with the Knight shift data displaying no decrease in spin 
susceptibility below T c for any field direction. 65 ' 66 The 
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combination of these features is incompatible with the 
usual pairing states such as the s-wave, p-wave, or d- 
wave state, and calls for an extension of the standard 
working scheme. 

In rcf. 29, we have investigated the magnetic prop- 
erties of non-centrosymmetric superconductors. Then, it 
was shown that the predominantly p-wave state admixed 
with the s-wave order parameter (s+P-wave state) is 
consistent with the paramagnetic properties of CePt 3 Si. 
We here examine the symmetry of superconductivity in 
CePt 3 Si from the microscopic point of view and show 
that the s+P-wave state or p+D+/-wave state can be 
stabilized by spin fluctuation with helical anisotropy. We 
also calculate the quasiparticle excitations, specific heat, 
and NMR 1/T\T, and show that the line node behav- 
ior in CcPt 3 Si at ambient pressure is consistent with the 
s+P-wave state. We investigate the pressure dependence 
of these quantities, possible SC multiple phase transi- 
tions, and the anisotropy of iJ C 2- Some future experi- 
mental tests are proposed. 

The paper is organized as follows. In §2, we formulate 
the RPA theory in the Hubbard model with ASOC and 
AFM order. The nature of spin fluctuation and super- 
conductivity is investigated in §3 and §4, respectively. 
The symmetry of superconductivity, SC gap structure, 
specific heat and NMR 1/TiT, multiple SC phase tran- 
sitions, and anisotropy of H C 2 are discussed in §4.1, §4.2, 
§4.3, §4.4 and §4.5, respectively. Some future experiments 
are proposed in §4. The summary and discussions are 
given in §5. A derivation of ASOC in the periodical An- 
derson model and Hubbard model is given in Appendix. 

2. Formulation 

2.1 Hubbard model with ASOC and AFM order 

For the following study of superconductivity in 
CePt 3 Si, we introduce the single-orbital Hubbard model 
including the AFM order and ASOC 

H = E £ ^H, s %s + a H 9(k) ■ S(k) - £ £q ■ S Q (k) 

k,s k k 

where S(k) = J2 s ,s' °ss'c\ c%. s ' and &Q( k ) = 
E 8 , s ' s c k.s' with ° ss ' bem S tnc vector represen- 

tation of the Pauli matrix. rii, s is the electron number 
at the site i with the spin s. We do not touch the heavy 
Fermion aspect, i.e., the hybridization of conduction elec- 
trons with Cc 4/-electrons forming strongly renormalized 
quasiparticlcs. However, we consider the Hubbard model 
as a valid effective model for describing low-energy quasi- 
particles in the Fermi liquid state. 67 

We consider a simple tetragonal lattice and assume the 
dispersion relation as 

e(k) = 2t\ (cos fc x + cos k y ) + 4^ cos fc x cos k y 
+2t 3 (cos 2k x + cos 2k y ) + [2<4 + 4t5(cos k x + cos fc y ) 
+4t 6 (cos 2fc x + cos 2k y )] cos fc z + 2t 7 cos 2fc z - fj,, (2) 
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determine the chemical potential /i so that 
the electron density per site is n. By choos- 
ing the parameters as (ti, t2, £3, <4, <5, tg, t-j, n) = 
(1,-0.15,-0.5,-0.3,-0.1,-0.09,-0.2,1.75), the dis- 
persion relation cq. (2) reproduces the /3-band of 
CePt 3 Si, which has been reported by band structure 
calculation without the AFM order. 68-70 The Fermi 
surface of this tight-binding model is depicted in Fig. 1 
of ref. 28. We assume that the superconductivity in 
CePt 3 Si is mainly induced by the /3-band because the 
/3-band has a substantial Ce 4/-electron character 68 and 
the largest density of states (DOS), namely 70% of the 
total DOS. 69 

The second term in eq. (1) describes the ASOC that 
arises from the lack of inversion symmetry and is char- 
acterized by the vector g(k). Time reversal symmetry 
is preserved, if the g- vector is odd in k, i.e., g(—k) — 
—g(k). In the case of CePt 3 Si as well as of CeRhSi 3 
and CeIrSi 3 , the g- vector has the Rashba type struc- 
ture. 71 The microscopic derivation of the ASOC in the 
/-electron systems is given in Appendix. The ASOC 
in the periodic Anderson model as well as that in the 
Hubbard model originate from the combination of the 
atomic L-S coupling in the /-orbital and the hybridiza- 
tion with conduction electrons. Although the detailed 
momentum dependence of the g-vector is complicated 
(see eq. (A. 25)) and is difficult to obtain by band struc- 
ture calculations, at least from a symmetry point of view, 
g(k) = (—v y (k),v x (k),0)/v delivers a reasonable approx- 
imation, where f x . y (fc) = de{k)/dk x>y is the quasiparti- 
cle velocity. We normalize g(k) by the average velocity 
v [v 2 — J2k v x(k) 2 + v y (k) 2 ) so that the coupling con- 
stant a has the dimension of energy. This form repro- 
duces the symmetry and periodicity of the Rashba-type 
g-vector within the Brillouin zone. We choose the cou- 
pling constant a = 0.3 in the main part of this paper so 
that the band splitting due to ASOC is consistent with 
the band structure calculations. 69 

The AFM order enters in our model through the stag- 
gered field hq without discussing its microscopic origin. 
The phase diagram under pressure implies that the AFM 
order mainly arises from localized Ce 4/-electrons that 
have a character different from that of SC quasiparti- 
cles. The T c of superconductivity is slightly affected by 
the AFM order which vanishes at P ~ 0.6GPa, 58 ~ 60 in 
contrast to the other Ce-based superconductors. 72 The 
experimentally determined AFM order corresponds to 
/iq = hqx pointing in the [100] direction with a wave 
vector Q = (0, 0,7r). 57 For the magnitude, we choose 
I /iq I «C W where W is the bandwidth since the observed 
AFM moment <~ 0.16/ib is considerably less than the full 
moment of the 5/2 manifold in the Ce ion. 57 

The undressed Green functions for U — are repre- 
sented by the matrix form G(k, iu n ) — (iw„l — P(fc)) -1 , 
where 

G(k,i. n )=( Cf^\ ^f^)\ (3) 
^ G 2 {k +1 iu n ) G 1 {k +1 iu n ) J 
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where the chemical potential /1 is included. We 
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and 



spin indices as 



H(k) 




e(k + ) 



(4) 



with e(k) = e(k)& ( ^ + ag(k)a and k + = k + Q. The 
normal and anomalous Green functions G l {k,'\uo n ) are 
the 2x2 matrix in spin space, where w„ = (2n + l)nT 
and T is the temperature. 

2.2 Eliashberg equation 

We turn to the SC instability that we assume to arise 
through electron-electron interaction incorporated in the 
effective on-site repulsion U. The linearized Eliashberg 
equation is obtained by the standard procedure: 



AA PjSl;S2 (fc) — 



k',q,S3,S4 



Vp,q,s 1 ,S2,S3,s i {k : , k )lpq,s 3 ,s 4 (k ), (5) 



(k)A mi (k"), (6) 



where q = p (q = 3—p) for i = j (i ^ j), k" = k+{i — l)Q 
and 

n 

xG^ S4 (-k + (p-l)Q,-iuj n ) (p=l,2).(7) 

Here, we adopt the so-called weak coupling theory of su- 
perconductivity and ignore self-energy corrections and 
the frequency dependence of effective interaction. 67,73 
This simplification strongly affects the resulting tran- 
sition temperature but hardly affects the symmetry of 
pairing. 67 We denote the order parameter for the super- 
conductivity as A PjSliS2 (A:) (p = 1,2, si and s 2 are the 
spin indices), where A 1>Sl;S2 (fc) and A 2iSliSi2 (fc) describe 
the Cooper pairing with the total momenta (0,0,0) and 
(0,0,7r), respectively. The former is the order parame- 
ter for ordinary Cooper pairs, while the latter is that for 
7r-singlct and 7r-triplct pairs. These 7r-pairs are admixed 
with usual Cooper pairs in the presence of the AFM or- 
der. 74 



The effective interaction V, 



p,q,s 1 ,s 2 ,s 3 ,s 4 



(k,k') originates 

from spin fluctuations that we describe within the RPA 73 
according to the diagrammatic expression shown in 
Fig. 1. In the RPA, the effective interaction is described 
by the generalized susceptibility whose matrix form is 
expressed as 

X2{q+) 



-X0) 



(0)/ 



*™® \ (8) 

xi°\u) ' (8) 



where q + = q + Q. Hereafter, we denote the element of 
the 4x4 matrix, such as A — x[°\q) and Xi{q) using the 



( A TTTT 


A 111l 


A 1111 


^TTU 






A UU 


^TUT 


Mill 


^ITTT 


A 11U 


A i\i\ 


Mill 




V ^UTT 


A llU 




Mill 


) 



(9) 



The matrix clement of the bare susceptibility Xi°\q) is 
expressed as 

Xi,]l,s2, s3,si^ = 

k.LO n 

+Gl 4M (k + q,iw n )G 2 s2<s3 (k + Q,iu; n )}, (10) 

X2,ll,s2,s3,s4:(q) = 

-^[G! 4ill (fc + q,iu n )G 2 s2tS3 {k,iuj n ) 

fc,W„ 

+Gl 4M (k + q,iw n )Gl 2<s3 (k + Q,iu; n )}, (11) 

and the matrix U is obtained as 

/ -U \ 

U 

U 

\ -U J 

According to Fig. 1, the effective interaction is obtained 
as 

Vl,l,al,s2,a3,s4(fc) k ) = V2,2,sl,s2,s3,s4(fc, k ) 

= -[U'xi{k' - k)U'] s3iSlM!s2 + U s i, S 2,s3,si, ( 13 ) 

(k,k') 

= -[U'X2(k' -k)U'] s3Ms2 , (14) 

where 

/ -u \ 
U 
U 
\ -U J 



u 



(12) 



V 



(15) 



x-xooc 



Fig. 1. Diagrammatic representation of the pairing interaction 

V P . 

tcraction U. 



The linearized Eliashberg equation (eqs. (5)-(7)) al- 
lows us to determine the form of the leading pairing insta- 
bility, which is attained for the temperature at which the 
largest eigenvalue A reaches unity. Numerical accuracy 
requires, however, a different but equivalent approach. 
We perform the calculation at a given temperature, in 
our case T = 0.02, which is much lower than the Fermi 
temperature, and determine the most stable pairing state 
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as the eigenfunction of the largest eigenvalue. 67 The typ- 
ical eigenvalue at T = 0.02 and U = 4 lies at around 
A = 0.3 ~ 0.6. This means that the T c for U = 4 is lower 
than T = 0.02. However, the absolute value of T c is not 
important for our purpose, which is focused on the roles 
of the spin fluctuation, ASOC, and AFM order. We be- 
lieve that the qualitative roles of these aspects can be 
captured in this simple calculation. On the other hand, 
the absolute value of T c is significantly affected by the 
mode coupling effect, vertex corrections, strong coupling 
effect, and multi-orbital effect, which are neglected in 
our calculation. We leave more sophisticated calculation 
based on the multi-orbital model and beyond the RPA 
for future discussion. 

The T c of superconductivity reaches T c — 0.02 if we 
assume a larger U. However, we show the results for 
U = 4, unless stated otherwise. This is mainly because 
the results for a large U are likely spurious because of 
the limitation of RPA. Since the mode coupling effect 
is neglected in the RPA, the critical fluctuation is not 
taken into account and therefore the magnetic instabil- 
ity is seriously overestimated. When we assume a large 
U so that we obtain a high T Cl the system approaches 
the magnetic instability, which is beyond the applicabil- 
ity of RPA. Our results for the superconductivity are 
only weakly dependent on U except for the relative sta- 
bility between the s+P-wave and p+.D+/-wave states, 
and therefore we qualitatively obtain the same results 
for a larger U. However, we avoid parameters close to 
the magnetic instability. 

Note that we here consider the spin fluctuation arising 
from quasiparticles that are mainly superconducting and 
may be different from the main source of the AFM mo- 
ment. Although the quantum critical point of the AFM 
order exists at P ~ 0.6GPa, the critical fluctuation of the 
AFM moment slightly affects the T c of superconductivity, 
as indicated by the phase diagram in the P-T plane. 58-60 
Therefore, it is expected that the critical fluctuations of 
the AFM moments are only weakly coupled to quasipar- 
ticles and the superconductivity is mainly induced by the 
residual interaction between quasiparticles. 

3. Spin Fluctuation 

First, we investigate the spin fluctuation in the non- 
centrosymmetric system. To clarify the role of ASOC, we 
consider the paramagnetic state where /iq = 0. Then, the 
static spin susceptibility is obtained as 

X ^(q) = [ dr < T T S"tf, r)S v {-q) > (16) 
Jo 

= CT sl,s2Xl,sl,a2, S 3, S 4(£)<3, S 4> ( 17 ) 

sl,s2,s3,s4 

where S^(q,r) = e HT S^(q) e - HT and S^(q) = 
Sfc, s ,s' a ss' c % + f s c k s '- We define Xma X (<?) as the maximum 
eigenvalue of the 3x3 matrix x^io) f° r M an d v. 

In the absence of ASOC, the spin susceptibility is 
isotropic, namely x^io) — for fi ^ v and x xx (o) = 
X vv (q) = X zz (q) = Xmax(<z)- The spin susceptibility has 
a peak at q = (0,0,7r) because of the band structure of 
the /3-band, as shown in Fig. 2(a) (dashed line). Thus, 



the /3-band favors the ferromagnetic spin correlation in 
the a&-plane and the AFM correlation between the plane. 
This is the spin structure realized in the AFM state of 
CePt 3 Si. 



(a) 

3 



o; 2 

x 

* E1 



--a=0 
— a=0.3 
-a=0.6 




Fig. 2. (color online) (a) Momentum dependence of spin suscep- 
tibility Xma,x{q) in the paramagnetic state (/iq = 0). We assume 
(7 = 4 and T = 0.02 and show the results for a = (dashed 
line), a = 0.3 (solid line) and a = 0.6 (dash-dotted line), (b) 
Spin susceptibility Xmax(9) along (0, 0, 7r)-(7r, 0, 7r) direction. We 
assume U = 3.95 for a = in (b) so that the maximum spin 
susceptibilities are similar between a = and a = 0.3. Numeri- 
cal calculation is carried out by dividing the first Brillouin zone 
into a 128 X 128 X 48 lattice. 



The anisotropy of spin susceptibility is induced by the 
ASOC. Our numerical calculation accurately takes into 
account the ASOC, but we explain here the role of ASOC 
within the first order of a to provide a simple and qual- 
itative understanding of the helical anisotropy. 

The lowest-order term in a appears in the off-diagonal 
component of the spin susceptibility tensor x^io) which 



is nonzero unless 



<7|| 



(n x ,n y )7r with n x and n y 



tegers. The off-diagonal component can be viewed as 
a result of the Dzyaloshinski-Moriya-type interaction 
Hum = T,q iD (v) ■ S(q) x S(-g). 41 Since the Rashba- 
type ASOC leads to D{q) oc az x q = a(—q yi q x ,0) in 
the vicinity of q^ = (0,0), the off-diagonal components 
of the spin susceptibility tensor are described as x yz (g) = 

-x zx (q) = 

0(a 2 ). Be- 



- X z Hq) - iaB(q z )q y + 0(a 3 ), X xz (g) = 
iaB(q z )q x + 0(a 3 ), and X xy (<l) = X yx (<7~) 



cause the momentum dependence of the diagonal com- 
ponent is quadratic as shown by X^^ (<lh 1z) ~ x(0,<? z ) — 
-<4(<Zz)|9|| | 2 + 0(|gj|| 4 ), the maximum eigenvalue of the 
spin susceptibility tensor is obtained as Xmax(<7y , Qz) — 
X(0,q z ) + aB(q z )\q {] \ - A(g z )|gj|| 2 + 0(a 2 , |g)|| 4 ) around 
qti = (0,0). Thus, Xmax(<7j|> <7z) for each q z has a local 
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minimum at gy = and a local maximum at gji ^ 0. For 
a = 0.3 the numerical calculation shows four peaks of 
Xmax(<zj|,<7z) at q ~ (0,±0.2,7r) and tf ~ (±0.2,0,tt) in 
contrast to the single peak at q = (0, 0, tt) for a = (see 
Fig. 2(b)). 

Since the off-diagonal components, such as x xz (?) = 
(x zx (<z))* an d X yz (<?) = (x zy (?))*! ar G purely imaginary, 
the maximum eigenvalue of the spin susceptibility tensor 
has the eigenvector S(q) = ^(g x ,g y ,±i) with q x _ y — 

9x,y/|<7)||- Thus, Xmax(<?|| i 9z) describes the susceptibility 
of the helical magnetic order. 

We here discuss the effect of higher-order terms of 
a by which the helical spin structure is distorted. In 
our result for a = 0.3, the spin susceptibility tensor at 
q = (0.196,0, 7r) has a maximum eigenvalue for S(q) ~ 
(0.81,0, -0.59i). The deviation from S(q) = ^=(1,0, ±i) 
in the lowest-order theory of a mainly arises from the 
second-order term of a in the diagonal component of 
X ^(q). 

Helical magnetism is suppressed by symmetric spin- 
orbit coupling, namely, the atomic L-S coupling, which 
is not taken into account in this paper. This is the reason 
why the helical magnetic order is not actually realized 
in non-centrosymmetric heavy fermion compounds, but 
is observed in non-centrosymmetric compounds with a 
small L-S coupling. 75 However, qualitatively the same 
effects of the ASOC, such as the helical anisotropy of spin 
susceptibility for gy ^ 0, are expected in the presence of 
L-S coupling. 

4. Superconductivity 

4-.1 Pairing symmetry 

We examine here the superconductivity. First, we dis- 
cuss the symmetry of the SC state. It is convenient in the 
following discussions to describe the order parameter in 
a standard manner as 76,77 



Ai )5 , s /(fc) 



-d x (k) 
-$(£?) 



idy(k) 
d z (k) 



$(fc) 
d x (k) 



d z (k) 
id y (fe) 



(18) 



where we use the even parity scalar function <I>(fc) and the 
odd parity d- vector d(k). In the presence of the AFM or- 
der, the order parameter for the 7r-triplet and 7r-singlct 
pairings A2, SjS '(^) appears owing to the folding of the 
Brillouin zone. However, the basic properties and sym- 
metries are hardly affected by 7r-parings when h,Q <C W. 

We identify two stable solutions of the Eliashberg 
equation. One pairing state has a predominant p-wave 
symmetry whose order parameter has the leading odd 
parity component d(k) ~ (— sin fc y , /3sin fc x , 0). The pa- 
rameter (3 is unity in the absence of the AFM order. The 
admixed even parity part is approximated as <$>{k) ~ 
S + cosfc x + cosfc y with S ~ 0.2. Thus, the spin sin- 
glet component has the s-wave symmetry, as discussed in 
ref. 24, but its sign changes in the radial direction in or- 
der to avoid the local repulsive interaction U . We denote 
this pairing state as the s+P-wave state. 

The other stable solution is the predominantly d-wave 
state that can be viewed as an interlayer Cooper pairing 



state: <&(&;) ~ {sin fc x sin fc z , sin fc y sin fc z } (two-fold degen- 
erate) admixed with an odd-parity component d(k) ~ 
<j>(fc)(— sin fc y , sin k x , 0). In the paramagnetic phase, the 
most stable combination of the two degenerate states is 
chiral: < &±(fc) ~ (sinfc x ± isinfc y )sinfc z which gains the 
maximal condensation energy in the weak-coupling ap- 
proach. In the AFM state, however, the two states of 
$(fc) are no longer degenerate. According the the RPA 
theory that we adopt in this paper, the d xz -wave state 
(d yz -wave state) is favored by the AFM order along the 
x-axis (y-axis) . Since the spin triplet order parameter has 
both the p-wave and /-wave components, we denote this 
state as the p+-D-l-/-wave state. 

In the RPA theory, superconductivity is assumed to 
be induced by the spin fluctuation. As discussed in §3, 
the spin fluctuation arising from the /3-band has four 
peaks around q ~ (0, ±0.2, tt) and q ~ (±0.2, 0, tt), which 
indicates the nearly ferromagnetic (helical) spin corre- 
lation in the a&-plane and the AFM coupling between 
the planes. For a small U, the spin fluctuation has a 
two-dimensional nature because of the dispersion rela- 
tion eq. (2). Then, the interplane AFM coupling is neg- 
ligible and the intraplane nearly ferromagnetic correla- 
tion induces s±P-wave superconductivity. On the other 
hand, the AFM coupling between the planes leads to a 
three-dimensional spin fluctuation for a large U and fa- 
vors the p+D+/-wave state. Figure 3 shows the phase 
diagram against U and the AFM staggered field /iq. We 
identify the pairing state with the largest eigenvalue in 
the Eliashberg equation at T = 0.02, as mentioned in §2. 
We see that the two pairing states are nearly degener- 
ate at around U = 3 ~ 3.5 independently of the AFM 
staggered field. 

We find the other pairing state having the predomi- 
nantly extended s-wave symmetry with <j>(fc) ~ cosfc x + 
cos fc y — 2 cos fc z as a self-consistent solution of the 
Eliashberg equation. However, we have found no param- 
eter set where this pairing state is stable. 



0, 



h 1 r 



p+D+f-wave 



s+P-wave 



0.05 0.1 0.15 0.2 

h Q 

Fig. 3. (color online) Phase diagram against U and /iq for a = 
0.3. The p+D-)-/-wave state (s+P-wave state) is stable for a 
large (small) U. The solid line shows the boundary of two phases. 
The Eliashberg equation is solved in the 128 X 128 X 32 lattice. 



Next we discuss the stability of the SC state when the 
ASOC is introduced. Figure 4 shows the a-dependence 
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of the eigenvalue of the Eliashberg equation A, for the 
s+P-wave and p+D+f '-wave states. The ASOC has two 
effects, namely, (i) the spin splitting of the band and (ii) 
the pairing interaction. The former is quantitatively im- 
portant in most non-centrosymmetric superconductors 
where \a\ > T c . It has been shown that the depairing 
effect due to (i) is minor for the spin singlet pairing 
state as well as for the spin triplet one with d(k) \\ g(k), 
while other spin triplet pairing states are destabilized. 2 
This is the reason why the s+P-wave state having the 
leading order parameter d(k) ~ (— sin fc y , /3sin fc x , 0) is 
favored among p-wave states that have a six-fold de- 
generacy in the absence of the ASOC and AFM order. 
Although the depairing effect due to the ASOC is al- 
most avoided in this s+P-wave state, it does not van- 
ish because the relation d(k) || g(k) is not strictly sat- 
isfied in the entire Brillouin zone. It should be noted 
that the momentum dependence of the vector in the ir- 
reducible represantation of the point group is not unique. 
Actually, the momentum dependences of the d-vector 
d{k) ~ (— sin fcy, sin k x , 0) and the g-vector g{k) = 
(— fy(fc), v K (k), 0)/v are inequivalent although these vec- 
tors are in the same irreducible representation of the Ci v 
point group. Although the momentum dependence of the 
d- vector is assumed to be the same as that of the g- vector 
in many theories, 21, 24-27,30-32,53 this assumption is not 
supported by the microscopic theory since the momen- 
tum dependence of the d- vector is mainly determined by 
the pairing interaction. Thus, Fig. 4 shows a steep de- 
crease in A in the s+P-wave state when a is turned on. 
This decrease arises from the depairing effect due to (i), 
whereby changes in the DOS due to band splitting is an 
additional source of the a-dependence. 




Fig. 4. (color online) Eigenvalues of Fliashberg equation A for 
the s+P-wave (circles) and p+D+/-wave states (triangles). We 
fix U = 4 and /iq = 0. The Fliashberg equation is solved in the 
128 X 128 X 32 lattice. The dashed lines show A which is estimated 
by using the spin susceptibility Xi(q) for a = 0. 



The effect (ii) of ASOC on the pairing interaction orig- 
inates from the modification of the spin fluctuation. This 
effect may be important in heavy Fermion systems since 
a large ASOC is likely induced through a strong L-S 
coupling in /-orbitals. Figure 2(a) shows the suppression 
of the spin susceptibility at around q = (0,0,7r), while 



that for other momenta is almost unchanged. Since the 
spin fluctuations around q — (0, 0, tt) are the main source 
of the pairing interaction in the p+D+f-w&ve state, 
the eigenvalue A for the p+P+/-wave superconductivity 
monotonically decreases asA = A(a = 0) — Aa 2 /e|. In 
contrast, the eigenvalue A for the s+P-wave state shows 
a minimum and increases with increasing a for a > 0.4. 
Thus, the effect of ASOC on the pairing interaction fa- 
vors the s+P-wave state rather than the p+D+f-w&ve 
state. 

These contrasting effects of ASOC arise from the 
anisotropy of helical spin fluctuation. To clarify this point 
we show the eigenvalue of Eliashberg equation A (dashed 
lines) in which the spin susceptibility Xi(q) is estimated 
for a — 0. We see that the A is increased for the s+P- 
wave state by the modification of spin susceptibility due 
to ASOC. Thus, the s+P-wave state is favored by the 
modified spin fluctuation although the magnitude of spin 
fluctuation is suppressed at around q = (0, 0, tt) by the 
ASOC. These results indicate that the anisotropy of he- 
lical spin fluctuation enhances the s+P-wave state. 

We confirmed that the subdominant component of the 
order parameter (s-wave component in the s+P-wave 
state, p- and /-wave components in the p+P^+/-wave 
state) grows almost linearly with increasing a, and that 
the momentum dependence of each component is almost 
independent of a. 

The eigenvalue of the Eliashberg equation is decreased 
by the AFM order owing to the loss of quasiparticle 
DOS. Therefore, the superconductivity is suppressed by 
the AFM order independently of the pairing symmetry. 
The relative stability of the s+P-wave and p+D+f -wave 
states is hardly affected by the AFM order, as shown 
in Fig. 3. The stability of the interlayer c?-wave state 
against the A-type AFM order has been claimed 78 by 
assuming the quasi- two-dimensional Fermi surface. How- 
ever, this is not the case in our model that assumes the 
three-dimensional /3-band. The T c of CePtaSi decreases 
if the AFM order decreases upon the application of pres- 
sure 58-60 and seems to be incompatible with our result. 
However, this pressure dependence may be due to the 
suppression of electron correlation by increasing pres- 
sure. 

It is expected that the AFM order leads to much more 
significant depairing effects on the intralayer d-wave and 
interlayer p-wave states because these Cooper pairings 
are directly broken by the A-type AFM order. The sta- 
bility of the SC state against the AFM order 58-60 implies 
the interlayer p+D+f -wave or intralayer s+P-wave state 
in CePtaSi that is identified in our calculation. 

4-2 Superconducting gap 

We investigate here the gap structure of both the s+P- 
wave and p+D+f-w&ve states and discuss the consis- 
tency with the line node behavior observed in CePtaSi 
at ambient pressure. 61-64 

The quasiparticle spectrum in the SC state is obtained 
by diagonalizing the 8x8 matrix using 



H s (k) = 



H(k) 
-AoAt(fc) 



-A A(fc) 



(19) 
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where A(k) is the SC order parameter in the spin basis 
expresses as 



(a) V° 



A(fc) 



Ai, S) y(fc)^ 



A 2iS) y(fc)^ 
A w (fc + Q) 



(20) 



In the following calculations the matrix clement of A(fc) 
is determined from the linearized Eliashberg equation by 
assuming that the momentum and spin dependences of 
the order parameter are weakly dependent on temper- 
ature for T < T c . We solve the Eliashberg equation at 
T = 0.02 > T c , having confirmed that the matrix A(k) 
is almost independent of temperature for T < 0.1. The 
same assumption has been adopted in other studies of 
multi-orbital superconductivity 79, 80 

Since the amplitude of the SC order parameter A is 
arbitrary in the linearized Eliashberg equation, we here 
choose Ao so that the magnitude of the maximal gap 
is A g = 0.1 in our energy units. Although this mag- 
nitude may be large compared with the energy scale a 
or /iq, we adopt this value for numerical accuracy, hav- 
ing confirmed that the lower values of A g do not alter 
the result qualitatively. We define the quasiparticle DOS 
p(e) as p(e) = j^J2 t J2'k S ( £ - E l {k)) where Yl k denotes 
the summation within the range \k z \ < ^. The eigenval- 
ues E 8 (k) > E 7 (k) > > E\{k) satisfy the relation 

Ei{k) = -Eg-i(k). 

It is more transparent to describe the SC order pa- 
rameter in the band basis, which is obtained by unitary 
transformation using 



A haDd (k) = U\k)A(k)U*(-k). 



(21) 



The unitary matrix U(k) diagonalizes the unperturbed 
Hamiltonian as 

U\k)H{k)U{k) = (e l {k)5 l] ). (22) 

The SC gap in the 7-th band is obtained as 

A 7 (fc) = A * 7 (fc), (23) 

where v E , 7 (fc) is the (77) component of the matrix 
Aband(^)- Since the relation T c <C |a| is satisfied in most 
of the non-centrosymmetric superconductors, the rela- 
tion |A b J and (fc)| < \ag(k)\,max{h Q ,\e(k) - e{k + Q)\} 
is valid for each (ij) component of Ah a nd(fc) except for 
the special momentum such as k = (0,0, k z ). There- 
fore, the off-diagonal components of Aband(fc) hardly af- 
fect the electronic state, and the quasiparticle excitations 
Ei(k) are approximated to ±E^ nd (k) with E^ nd (k) 2 = 

e 1 (k) 2 + |A 7 (fc)| 2 . Thus, the SC gap in the 7-th band is 
described by |A 7 (fc)|. 

It is clear that the p+D-|-/-wave state has a horizon- 
tal line node protected by the symmetry because all of 
the matrix elements of A(fc) are zero at k z = 0. This 
is consistent with the experiments in CePtaSi. 61-64 The 
coefficient c\ of the linear term in the DOS (p(s) = cie) 
increases in the AFM state because the pairing state 
changes from the chiral d xz ± ic? yz -wavc state in the para- 
magnetic state to the d xz -wave state in the AFM state 
(see Fig. 3 of ref. 28). 




(b) 



oc=0.3 



0.5 




Fig. 5. (color online) DOS p(e) in the s+P-wave state at U = 4. 
(a) Paramagnetic state for a = (dashed line) and 0.3 (solid 
line), (b) AFM state with ASOC a = 0.3 for h Q = 0.125 (solid 
line) and h,Q = 0.2 (dash-dotted line). We show the results for 
e > because p(e) is particle-hole symmetric owing to its defini- 
tion. The Eliashberg equation is solved in the 64 X 64 X 32 lattice 
and the DOS is calculated in the 384 X 384 X 384 lattice. 



We investigate here the accidental line node of the SC 
gap in the s+P-wave state in detail. The quasiparticle 
DOS p(e) is shown in Fig. 5 and the SC gaps |A 7 (/c)| 
for 7 = 3, 4 are shown in Figs. 6 and 7, respectively. The 
paramagnetic state is assumed in Fig. 6, while the AFM 
state is assumed in Fig. 7. 

The SC gap in the absence of the ASOC and AFM 
order (Fig. 6(a)) is approximated to be |A(/c)| ~ 



sin fc 2 + sinA; 2 , which has two point nodes in the [001] 

direction in contrast to the experimental results. 61-64 
The DOS at low energies is quadratic as shown by 
p(e) ~ C2S 2 and the coefficient c 2 is small owing to 
the small DOS in the [001] direction (dashed line in 
Fig. 5(a)). 

The line nodes are induced by the ASOC through the 
following two mechanisms. 

(I) Admixture with an s-wave order parameter. 

(II) Mismatch of the d-vector and g-vector. 

The first one (I) has been proposed by Frigeri et al. A1 
and its contributions to the NMR \jT\T and superfluid 
density have been investigated by Hayashi et al. 30 ' 32 In 
the absence of the AFM order, the SC gap is expressed 
as ±$(fc) + d(k) ■ g(k) with g(k) = g{k)/\g(k)\. The s- 
wave and p-wave order parameters are approximated to 
be <&(fc) ~ <&(0, k z ) and d(k) ■ g(k) ~ c(fc z )|fc||| at around 
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Fig. 6. (color online) SC gap in the s+P-wave state at k z = j in 
the paramagnetic state (/iq = 0). (a) SC gap in the absence of 
the ASOC (a = 0). We show the SC gaps |A 3 (fc)| and |A 4 (fc)| 
for a = 0.3 in (b) and (c), respectively. The SC gaps [Ai(fc)| and 
|A2(fc)| arc not shown because the quasiparticlc DOS is small 
in these bands. The positions of the line nodes arising from the 
mechanisms (I) and (II) are shown by the arrows in (b) (see the 
text). The thin dashed lines show the Fermi surface at fc z = 
Although the zeros of SC gap do not intersect with the Fermi 
surface at k z = ~, the line nodes exist on the Fermi surface 
at another k z . Note that the /3-band has a three-dimensional 
Fermi surface, while the SC gap in the s+P-wave state is nearly 
independent of fc z . 



k\\ = (k K ,k y ) = (0,0), respectively. Therefore, the SC 

gap vanishes on the line |fc||| = |$(0, fc z )|/c(fc z ) in half 
of the bands, while the other bands have a full gap. We 
show the SC gaps |A3(fc)| and |A4(fc)| at fc z — tt/3 in 
Figs. 6(b) and 6(c), respectively. The line node actually 
appears in |A 3 (fc)| in the vicinity of fc|| = (0,0) (shown 
by the arrow (I)). However, the line node arising from the 
mechanism (I) induces only a tiny linear term p(e) ~ c±e 
with ci oc | a | because the length of the line node is very 
small, as shown in Fig. 6(b). 

We find another line node arising from the mecha- 
nism (II) at around |fcn| = 7r/3 (see the arrow (II) in 
Fig. 6(b)). This line node originates from the topological 
character of the g- vector. According to the assumptions 
9x(k) = —v y (k)/v and g y {k) — v K (k)/v, the g-vector 
has a singularity not only on the [001] line but also on 
the line at around (k x ,k y ) = (0.47T, OAir). The g-vector 
rotates around the singular point, and therefore the re- 
lation d(k) _L g(k) is satisfied on a line. The SC gap 
±$(fc) +d(k) ■ g(k) vanishes around this line because the 
s-wave component |3>(fc)| is much smaller than the p- 
wave component |d(fc)|. This is a general mechanism for 
the line node in the non-centrosymmetric superconduc- 
tor predominated by the spin triplet pairing. However, 
it is not clear whether this line node exists in CePtaSi 
because it depends on the detailed momentum depen- 
dence of the g- vector. For example, this line node does 
not appear if we assume g(k) = (— sin k y , sin fe x , 0). Any- 
way, the low-energy excitation arising from the ASOC 
is small because of the steep increase in SC gap around 
the line node, as shown in the schematic figure (Fig. 8). 
Figure 5(a) actually shows a small coefficient c\ of the lin- 
ear term p(e) — c±e (solid line). This linear term mainly 
arises from the line node induced by the mechanism (II) 
and the contribution of the line node (I) is negligible. 

The DOS at low energies is markedly increased by the 
AFM order owing to the following two effects. 

(III) Folding of the Brillouin zone. 

(IV) Mixing of the p-wave order parameter between the 
leading part d{k) <~ (— sin k y , sin fc x , 0) and the admixed 
part d(k) = (sin k y , sin fc x , 0). 

The former arises from the pair-breaking effect due to 
the band mixing, which has been investigated by Fuji- 
moto. 31 In contrast to ref. 31, the line node appears not 
only at k z — n/2 but also at around (fcx, k y ) — (n/2, it/2) 
(see Fig. 7(b)) in our case because of the band structure 
of the /3-band. However, the DOS arising from (III) is 
not quantitatively important when /iq <C W because of 
the steep increase in the SC gap around the line node, 
as shown in Fig. 8. 

Actually, the low-energy excitations in the s+P-wave 
state are mainly induced by the effect (IV). The a- and b- 
axes in the tetragonal lattice are no longer equivalent in 
the presence of the AFM order. Therefore, the p-wave or- 
der parameter is modified to d{k) = (— sinfc y ,/3sinfc x ,0) 
with (3 =/= I. This change can be viewed as the mixing 
of the leading part d(k) = (— sin fc y , sin fc x , 0) with the 
admixed part d(k) = (sin k y , sin fc x , 0), which leads to 
the rotation of the d-vector. According to the result ob- 
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Fig. 7. (color online) SC gap in the s+P-wave state with the 
AFM order. We show (a) |A 3 (fe)| and (b) |A 4 (fe)| at k z = f for 
a = 0.3 and h,Q = 0.125. The positions of the line nodes arising 
from the mechanism (III) and (IV) are shown by the arrows in 
(b). 



tained using the RPA theory, f3 decreases with increasing 
ft.Q. Then, many low-energy excitations are induced at 
around k y — n/6, as shown in Fig. 7(b). The SC gap in 
the 4-th band (Fig. 7(b)) is further decreased at around 
fc y = 7r/6 by the admixture with an s-wave order pa- 
rameter. The DOS clearly shows a linear dependence in 
Fig. 5(b), which is consistent with the experimental re- 
sults in CePtaSi at ambient pressure. 6 1-64 We have shown 
that the rotation of the d- vector is also the main source 
of the anomalous paramagnetic properties of CePtaSi. 29 

4.3 Specific heat and NMR 1/T X T 

The pressure dependence of the SC state is a deci- 
sive test for validating the theory of CePtaSi as well as 
of CeRhSi3 and CeIrSi3. According to the experimental 
result of CePt 3 Si, 58-60 the AFM order is suppressed at 
a pressure P ~ 0.6GPa, although the superconductiv- 
ity survives at high pressures P > 0.6GPa. Therefore, 
the role of the AFM order can be studied experimentally 
by measuring the pressure dependence of the SC state. 




Momentum 



Fig. 8. (color online) Schematic figure showing the momentum 
dependence of SC gap along the Fermi surface. The solid line 
shows the usual line node protected by the symmetry (for exam- 
ple, the d x 2_ y 2-wave superconductor). The accidental line nodes 
(II) in Fig. 6(b) and (III) in Fig. 7(b) show a steep increase in 
the SC gap around the gap node (dashed line). 



If the s+P-wave state is realized in CePtaSi and the 
AFM order is the main source of line nodes, the number 
of low-energy excitations decreases under pressure. This 
theoretical result can be tested by measuring the pres- 
sure dependence of specific heat, NMR l/T^T, superfluid 
density, thermal conductivity, and other quantities. We 
now calculate specific heat and NMR \jT\T for a future 
experimental test. 

To discuss these quantities, we adopt the same assump- 
tion in §4.2. We here calculate the amplitude of the SC 
gap, A n , in eq. (19) by solving the gap equation 

1 = gJ2 l* 7 (fc)| 2 tanh -^^/2E^ d (k), (24) 
fc 

which is obtained as a mean field solution of the effective 
model in the band basis given as 

H = E E ^) d \/kn - ^E b l b *> (25) 

k 7—1 k,k' 
7 

The SC order parameter obtained in the linearized 
Eliashberg equation (eqs. (5)-(7)) is reproduced using 
this model. We choose g so as to obtain T c = 0.05. We 
have confirmed that the smaller g and T c do not quali- 
tatively alter the following results. 

The quasiparticle excitation Ei(k) is determined using 
eq. (19) with A determined using eq. (24). The Som- 
merfeld coefficient C/T is obtained as 

C / T =qT> (27) 

s = - EE[/a lQ s fit + ^ - 4s) ^(i - 4s)]. 

k i=l 

(28) 

where /. g is the Fermi distribution function f-t— (1 + 

exp(S 4 (fc)/T))- 1 . 

We calculate NMR 1/TiT as 

i/r 1 r = imxL(n)/n| t wo, (29) 
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XL(i^n) = - ^2 ^2\Gi^{k' ,iuj n + in n )G l:ll (k,iuj n ) 

k,k',uj n i—1 

-^ iT (fe',iw„ + in n )F UT (fc,icj n )], (30) 

where Gi(k,iuj n ) and Pj(fc,iaj„) are the normal and 
anomalous Green functions in the SC state, respectively. 
We ignore the momentum dependence of the hyperfine 
coupling constant and the exchange enhancement due to 
the electron correlation for simplicity. The local spin sus- 
ceptibility xl(^) is obtained from Xh(^n) through the 
analytic continuation. 



NMR \jT\T shows a T 2 dependence at low tempera- 
tures. These results are consistent with the experimental 
data of specific heat, 63 thermal conductivity, 61 superfluid 
density, 62 and NMR l/T^T. 64 

Upon decreasing the staggered field /iq, low-energy 
excitations are suppressed. In the paramagnetic state 
(Hq = 0), the Sommerfcld coefficient deviates from the 
T-linear law below T < 0.2T C , while the T 2 depen- 
dence of NMR l/TiT breaks down at lower tempera- 
tures, T < 0.1T C (solid lines in Fig. 9). If the AFM order 
is the main source of the line node in CePtaSi at ambi- 
ent pressure, these deviations from the line node behavior 
may be observed at high pressures P > 0.6GPa. 
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Fig. 9. (color online) (a) Sommerfeld coefficient C/T and (b) 
NMR 1/TiT in the s+P-wave state. The parameters a and /iq 
are shown in the figure. We solve the Eliashbcrg equation in 64 
X 64 X 32 lattices and estimate C/T and 1/TiT in 384 X 384 
X 192 lattices. 



Figures 9(a) and 9(b) respectively show the tempera- 
ture dependences of the Sommerfeld coefficient C /T and 
the NMR 1/TiT in the s+P- wave state. When we as- 
sume a weak ASOC (\a\ <C T c ) and the absence of the 
AFM order (dashed lines in Fig. 9), both the Sommerfeld 
coefficient and the NMR \jT\T at low temperatures are 
much smaller than those expected in the superconductor 
with line nodes. For example, the Sommerfeld coefficient 
shows a T" dependence (n > 2) which is incompatible 
with the experimental result. 63 On the other hand, we 
clearly see the line node behavior in the presence of the 
ASOC and AFM order (dash-dotted lines in Fig. 9). The 
Sommerfeld coefficient obeys the T-linear law and the 
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Fig. 10. (color online) (a) Sommerfeld coefficient C/T and (b) 
NMR 1/TiT in the p+D+/-wave state. 



Figure 10 shows the Sommerfeld coefficient and NMR 
1/TxT in the p+D+f -wave state. The line node behav- 
ior appears clearly in both the paramagnetic and AFM 
states. The role of the AFM order is qualitatively the 
same as that in the s+P-wave state: the number of low- 
energy excitations is increased by the AFM order. This 
is because the vertical line node in the d xz -wave state 
disappears in the chiral <i-wave state. 

We here discuss the coherence peak in the NMR 
\jT\T. It has been shown that the coherence peak ap- 
pears in the s+P-wave state just below T c owing to the 
finite coherence factor. 20 ' 30 This is the case in our cal- 
culation; however, the coherence peak is much smaller 
than that shown in ref. 30, as shown in Fig. 9(b). This 
is because of the small ASOC a = 0.3 <C £f assumed in 
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this paper and the extended s-wave nature of the spin 
singlet order parameter. The coherence factor in the ex- 
tended s-wave state is decreased by the sign reversal of 
the order parameter in the radial direction. Note that 
the isotropic s-wave pairing is generally not favored in 
the strongly correlated electron systems. A slightly larger 
coherence peak appears in the paramagnetic state (solid 
line in Fig. 9(b)); however, this is not due to the co- 
herence factor but arises from the anomaly in the DOS. 
Although a coherence peak was reported in the early 
measurement of NMR 1/TiT, 56 the recent measurement 
for a clean sample shows no coherence peak just below 
T c , 64 in agreement with our result. 

4-4 Multiple phase transitions 

We have discussed the pressure dependence of low- 
energy excitations in §4.2 and §4.3. Although qualita- 
tively the same results are obtained for the low-energy 
excitations between the s+P-wave and p+Z?+/-wave 
states, there is an essential difference, namely, the multi- 
ple phase transitions in the P-T plane. To illustrate this 
issue, we show the possible phase diagrams in Fig. 11. 



(a) p-wave (a » Tc) 
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(b) p-wave (a < Tc) 
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Fig. 11. (color online) Possible phase diagrams in the P-T plane, 
(a) s+P-wave state for a large ASOC {\a\ > T c ). (b) .s+P-wave 
state for a small ASOC (|a| < T c ). (c) p+D+/-wave state for 
a large staggered field Iiq. (d) p+P>+/-wave state for a small 
staggered field. "D" ("cD") shows the d xz -wave (chiral d xz ± 
idy Z -wave) state. "cP" shows the chiral p-wave state where the 
dominant order parameter is d = (p x ± ip y )x. 



Figures 11(a) and 11(b) show the phase diagrams in 
the s+P-wave state. When the ASOC is small (|a| < T c ), 
the chiral p-wave state is stabilized at low temperatures 
and low pressures, as in Fig. 11(b). However, this is un- 
likely for CePtaSi since the ASOC is much larger than 
T c in heavy fermion systems. Therefore, the simple phase 
diagram in Fig. 11(a) is expected in the s+P-wave state 
of CePt 3 Si. 

In the case of the p+D+f -wave state, the phase transi- 
tion from the chiral d xz ± ic? yz -wave state to the ci xz -wave 
state must occur, as in Fig. 11(c) or 11(d). When the 
staggered field /iq is large (small) at ambient pressure, 



the phase diagram in Fig. 11(c) (Fig. 11(d)) is expected. 
Thus, the enhancement of the low-energy DOS due to 
pressure accompanies the second order phase transition, 
in contrast to the s+P-wave state. The observation of a 
multiple phase transition in the P-T plane might provide 
clear evidence of the p+D+/-wave case. Although the 
second SC transition has been observed in CePtsSi, 62 ' 81 
it has been shown that there are two SC phases with 
T c - 0.75K and T c ~ 0.45K in the sample. 63 ' 81 ~ 83 The 
second transition below T c seems to be caused by sample 
inhomogeneity. 

4-5 Anisotropy of upper critical field H C 2 

We here comment on the in-plane anisotropy of H C 2 
arising from the AFM order. As discussed in J4.2, 
the p-wave order parameter in the s+P-wave d(k) ~ 
(— sinfc y ,/3sinfc x ,0) has a two-fold in-plane anisotropy 
in the AFM state. The anisotropy parameter (5 can be 
measured by the in-plane anisotropy of H C 2 near T c , 
which is determined by the orbital depairing effect and 
written as H* 2 = $ /(27r£ b £ c ) and H* 2 = $ /(27r£ a &) 
for H \\ a and H || b, respectively. Here, 3>o = ff 
is the flux quantum and ^ h , c = £° b>c (l - T/T,.)- 1 / 2 
are coherence lengths. We obtain the ratio of the gra- 



a, b ,c 



T c dH^' c /dT as, H c l : H c b 2 : H, 



'b 



c2 



dient H ( 

£a • £b ■ £c which can be estimated using the relation, 
(C°b,c) 2 « EU<, ai b,c(^)l^(fc)l 2 /"(e 7 (fc))/8 e7 (fc), 
where i> 7 , a , b ,c(fc) = de 7 (fc)/dfc aibjC is the quasiparticle ve- 
locity in the 7-th band. 

Figure 12 shows the in-plane anisotropy H c2 /H C 2 in 
the s+P-wave state (solid line). It is clearly shown that 
the anisotropy is induced by the AFM order for hq > 
0.1. This is mainly due to the decrease in the anisotropy 
parameter (3. Since (3 < 1 in the RPA theory and we 
assume the AFM staggered moment pointing along the a- 
axis, H C 2 is higher along the 6-axis than along the a-axis 
(H* 2 < H^ 2 )- If /? > 1, the opposite anisotropy appears. 
Thus, if the marked mixing of p-wave order parameters 
due to the AFM order occurs, a pronounced in- plane 
anisotropy appears in H c2 . 

The paramagnetic depairing effect qualitatively in- 
duces the same in-plane anisotropy as that in Fig. 12. 
We have shown in ref. 29 a schematic figure of the H-T 
phase diagram by taking into account both the orbital 
and paramagnetic depairing effects. 

The in-plane anisotropy of H c2 in the p+D+J-wave 
state is quite different from that in the s+P-wave state. 
We obtain HJ^/Hjz ~ 1.6 in the p+£)+/-wave state 
independent of the nonzero staggered field /iq. This 
two-fold anisotropy changes discontinuously way from 
H c2.l H 'd2 ~ 1-6 in the AFM state to H'£/H'£ = 1 in 
the paramagnetic state, in contrast to the continuous 
change in the s+P-wave state. It is therefore expected 
that the s+P-wave state can be distinguished from the 
p+-D+/-wave state by the pressure dependence of in- 
plane anisotropy in H c2 . 

Next we comment on the experimental measurement 
of in-plane anisotropy arising from the AFM order. The 
direction of the AFM moment can be controlled by the 
cooling process, namely zero- field cooling and field cool- 
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Fig. 12. (color online) In-plane anisotropy of H c 2 in the s+P- 
wave (solid line) and p+D+/-wave (dashed line) states. We 
assume a = 0.3. We show the ratio of the gradient H c2 = 
-T c dH* 2 /dT and H'J% = -T c dH^ 2 /dT where H* 2 and H^ 2 are 
the upper critical fields along a- and b-axes, respectively. 



ing. When temperature is decreased under the magnetic 
field along the 6-axis, the AFM moment parallel to the 
a-axis appears below the Neel temperature, because the 
system gains the maximum magnetic energy when the 
AFM moment is perpendicular to the magnetic field. 
Then, the two-fold anisotropy due to the AFM order ap- 
pears at low magnetic fields, although the AFM moment 
may rotate at high magnetic fields. On the other hand, 
the domain structure with respect to the direction of 
AFM moment can appear when the system is cooled un- 
der a zero magnetic field. Then the two- fold anisotropy 
is obscured. 

Before closing this section, some comments are given 
on the anisotropy of H C 2 between the a&-plane and the 
c-axis. We cannot discuss this anisotropy in a final way 
because not only the /3-band but also the other band 
affects the anisotropy. However, it should be noted that 
H C 2 is of similar magnitude along the a6-plane and the 
c-axis because the /3-band has a three-dimensional Fermi 
surface. For example, we obtain Pf c | : Hj£ : = 1:1: 
0.71 for hq = and H' c % : H% : H' c c 2 = 0.67 : 1 : 0.55 
for /iq = 0.125 in the s+P-wave state, while H c ?, : : 
H' c l = 1 : 1 : 0.85 for h Q = and H',% : H'£ : H' c c 2 = 
1.61 : 1 : 0.86 for h Q = 0.125 in the p+D+f-wave state. 
The weak anisotropy of H C 2 between H || ab and H \\ c 
is consistent with the experimental result for CePtaSi. 58 

5. Summary and Discussion 

We have investigated the superconductivity in the 
Hubbard model with Rashba-type spin-orbit coupling 
and AFM order. Applying the RPA theory to the /3- 
band of CePtaSi, we found two stable pairing states, the 
intraplane p-wave state admixed with the s-wave compo- 
nent (s+P-wave state) and the interplane d-wave state 
admixed with the p- and /-wave components {p+D+f- 
wave state) . We found that the anisotropy of helical spin 
fluctuation favors the s+P-wave state. 

We examined the low-energy excitations in detail. The 
SC gap in the p+D+f-wave state has a line node pro- 
tected by the symmetry, while accidental line nodes ap- 
pear in the s+P-wave state. Thus, both pairing states 
seem to be consistent with the experimental results in 
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CePtaSi at ambient pressure. 61-64 A substantial part of 
the accidental line node in the s+P-wave state can be 
induced by the AFM order through the rotation of the 
d-vector. The line node in the p+D+f-wave state is also 
increased by the AFM order because of the phase tran- 
sition from the chiral d-wave state in the paramagnetic 
state to the <i xz -wave state in the AFM state. Thus, the 
number of low-energy excitations decreases in both states 
when the AFM order is suppressed by pressure. We cal- 
culated the specific heat and NMR \jT{F in both the 
paramagnetic and AFM states. The deviation from the 
line node behavior in the paramagnetic state has been 
pointed out. 

We proposed some future experiments that can eluci- 
date the pairing state in CePtaSi. The first one is the 
pressure dependence of low-energy excitations discussed 
above. Another one is the possible multiple SC phase 
transitions in the P-T-plane. The second SC transition 
occurs below T c near the critical pressure for the AFM 
order, if the p+D+/-wave superconductivity is realized. 
This is in contrast to the s+P-wave state where no addi- 
tional phase transition is expected. The marked change 
of low-energy excitations in the p+D+f-wave state is 
accompanied by the second order phase transition. The 
last proposal is the anisotropy of H C 2 in the a&-plane. In 
the s+P-wave state, the anisotropy of H c2 gradually in- 
creases with increasing AFM moment, while that in the 
p+D+f-wave state is discontinuous at a critical pressure 
for the AFM order. Our proposals for future experiments 
do not rely on the particular band structure of the (3- 
band in CePtaSi, and therefore can also be applied to 
CeRhSia, CeIrSi 3 , and CeCoGe 3 . 

According to the present experiments, the s+P-wave 
superconductivity is most likely realized in CePtaSi. 
The paramagnetic properties measured on the basis of 
the NMR Knight shift and H c2 seem to be compatible 
with those in the s+P-wave state. 29 Futher studies from 
both the theoretical and experimental points of view are 
highly desired to elucidate the novel physics in the non- 
centrosymmetric superconductivity. 
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Appendix: Derivation of Rashba-type ASOC in 
Tight-binding Models 

We here microscopically derive the Rashba-type spin- 
orbit coupling in the periodical Anderson model and 
Hubbard model by tight-binding approximation. 

The localized 4/ states in the Ce-based heavy fermion 
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superconductors, such as CePtaSi, CeRhSi3, and CeIrSi3 
are described by the J = 5/2 manifold whose degener- 
acy is split by the crystal electric field. The 4/ levels in 
CePt 3 Si are described by the three doublets, 84 



ir 7 ± >= 



|r' 6 ± >= | ± - >, 




(A-l) 
(A-2) 

(A-3) 



The ground state is |T7± >, and the excited \T' 6 ± > and 
|r' 7 ± > states have excitation energies of 1 and 24mcV, 
respectively. 57 

Next we construct a periodical Anderson model for the 
|r7± > state, which hybridizes with conduction elec- 
trons. It is straightforward to apply the following pro- 
cedure to the |Tg± > and |T 7 ± > states. Because the 
mirror symmetry is broken along the z-axis in CePtsSi, 
the odd parity Ce 4 /-orbital is hybridized with the even 
parity s- and d-orbitals in the same Cc site. Owing to 
the symmetry of the |r7± > state, the localized 4/ state 
is hybridized with the d xy -, d KZ - and d yz -orbitals. Then, 
the wave function of the localized state can be expressed 
as, 



|f± >= K|r 7 ± > +ie|d xy > X ± 
+??(|d xz > Ti|d yz >)x T , 



(A-4) 



where e, rj and k = \J\ — e 2 — 2r\ 2 are real and \± de- 
scribes the wave function of the spin. We note that the 
wave function of |T7± > is given by 



|r 7 ± >= 



2- >X± 



±( V^I |jL " = ±3 > " V 2T |iz = Tl >)Xt ' 



(A-5) 



where \L Z = 2- >= ^(|L Z = 2 > -\L Z = -2 >). 

The periodical Anderson Hamiltonian is constructed 
for the localized |f± > state and conduction electrons. 
We here consider the conduction electrons arising from 
the Ce 5s-orbital for simplicity. Taking into account the 
inter-site hybridization between the s-, d- and /-orbitals, 
we obtain the tight-binding Hamiltonian 



H = E^o(%, 



where tbt=(fl , ft , ct , ct , ) , and 



(A-6) 



H {k) 



ii(k) 
t>(fc)t 



V(k) 
s c (k) 



(A-7) 




V(k) = 

(8V 3 s z + 4iet4)s x s y (2iV 5 - 4r]V 6 s z )(s X + lSyj 
(2iV 5 - AtjVgS^ (s x - is y ) (8V 3 s z + 4ieV" 4 )s x s y 

(A-10) 

where the abbreviation s Xiy , z = sin fc x ^ y ^ z is used. We ig- 
nored the off-diagonal terms in the second order with 
respect to the small parameters e and r\. We obtain 
£ f (fc) = K 2 e rr (k) + e 2 e xy (fc) + r) 2 (e xz (k) + e yz (fc)) where 
£A(fc) is the dispersion relation for the |A > state. It is 
clearly shown that eq. (A. 8) has the Rashba type spin- 
orbit coupling term and the coefficient is obtained as 



a 1 = -4eV 2 - ArjVi. 



(All) 



The hybridization parameters in cqs. (A. 10) and (A. 11) 
are obtained as 



Vi 

v 6 



— V lm 




T/100 

ir3 — 3x 2 y,yz 



ylOO 



2 ),yz 



(A-12) 



),(A13) 



V 



no 



xy, s; 



k(J 15 T/!00 



3xy 2 ,s 



— y 100 

42 x ( 5z 



2 ),s 



v 101 

v XZ.S 5 



(A-14) 
(A-15) 

. (A-16) 
(A-17) 



where is the hopping matrix element between the 

| A > and |B > states along the [abc]-axis. 

Note that the parameters e and rj arise from the intra- 
site hybridization between the d- and /-orbitals while 
the matrix elements V\ and V 2 describe the inter-site 
hybridization. Thus, the intra-orbital Rashba-type spin 
orbit coupling ol\ arises from the hybridization of the 
r7-state with the <i xy -, d xz -, and d yz -states. Note again 
that the parameters e and r\ vanish in centrosymmetric 
systems. 

Applying an appropriate unitary transformation to the 
conduction electron, (ct , ct ) = (ct , ct )U c (k), the hy- 
bridization matrix is transformed as 



V(k) - V(k)U c (k) 
V c{ 0) 



a 2 (fc)(is y s x + s 2 s y ) 
V ci {k) 



\ a 2 (fc)(-is y s x + s 2 s y ) 
where 

a 2 {k) = 4{eViV 5 - A V V 3 V 6 s 2 z ) 

I yj 'l6y 3 2 s2 S 2 S 2 + V£ (S x + S y ). 



,(A18) 



(A-19) 



The 2x2 matrix if(k), i c {k), and V(k) are obtained as 
ai(-is x + Sy) e f (fc) 



(A-9) 



Note that a 2 (k) is a real and even function with respect 
to A) x ; ky ; ciiid k z . 

Taking into account the on-site repulsion in the |f± > 
state, we obtain the periodical Anderson model with a 
Rashba-type ASOC as 



H = -ff k + -ffAsoc + Hi, 



(A-20) 



14 J. Phys. Soc. Jpn. 

fc,S=± fc,S=± 

+ E [Kf(fc)4 sCfe - s + fc. C .], (A-21) 

fc,s=± 

#ASOC = <*i E 5f(fc) ' °ss'fljk,s> 

k.s,s' 

+ E Mfc)&f(fc) • Z»'fl a %., + h-c], (A-22) 

fc,s,s' 

ffi = ^E<+<- ( A ' 23 ) 

z 

where (/f(fc) = (sin fc y , — sin fc x , 0) and g c f(k) = 
(sin 2 fc x sin fc y , — sin 2 fc x sin fc x , 0) describe the g- vector for 
the intra- and inter-orbital Rashba-type ASOCs, respec- 
tively. 

Note again that the ASOC arises from the atomic L-S 
coupling in the Ce 4/-orbital and the parity mixing in 
the localized state. The breakdown of the inversion sym- 
metry plays an essential role in the parity mixing in the 
atomic state. The inter-site hybridization between the 
/- and admixed c?-(or s-)orbitals gives rise to the intra- 
orbital ASOC, while the inter-orbital ASOC is induced 
by the hybridization between the conduction electrons 
and the admixed d- (or s-) orbit als. Note that the cubic 
term cx sin fc x sin k y sin k z (sin 2 fc x — sin 2 fc y )cr z 26 does not 
appear in the above derivation. 

We derived the periodical Anderson model for the 
localized |Tg± >= | ± | > state and the conduc- 
tion electrons with s-orbital symmetry. Then, we ob- 
tained the Hamiltonian that is similar to eq. (A. 20), 
but the (/-vector for the inter-orbital ASOC is replaced 
with g c f(k) — (sin fc y , — sin fc x , 0) = !Jf(k)- Other crys- 
tal field levels different from eqs. (A.1)-(A.3) have been 
proposed. 2 ' 85 The Rashba-type ASOC can also be de- 
rived for these levels in the same way as above. Thus, 
the Rashba-type ASOC is generally derived in the peri- 
odical Anderson Hamiltonian by taking into account the 
parity mixing in the atomic 4/-state. 

The kinetic energy term in the periodical Ander- 
son model is diagonalized by the unitary transformation 

^(£)=f 0l f§ (A24) 

Applying this unitary transformation to the periodical 
Anderson model in eq. (A. 20) and dropping the upper 
band described by ^ , we obtain the single-orbital 
model with the Rashba-type ASOC. The (/-vector is ob- 
tained as 

ag(k) = aiai(fc) 2 (7f(fc) + 

a 2 (k) ai (k)(a 2 (k) + a* 2 (k))g ct (k). (A-25) 

The unitary transformation described by U c f(k) leads to 
the momentum dependence of the two-body interaction 
term Hi, as in the case of the multi-orbital Hubbard 
model. 86 By neglecting this momentum dependence for 
simplicity, we obtain the single-orbital Hubbard model 
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in eq. (1) where the (/-vector is described by eq. (A. 25). 
The investigation of the periodical Anderson model in 
eq. (A. 20) is an interesting future issue. 
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